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This paper presents a useful compact formula for deriving an effective Hamiltonian describing 
the time-averaged dynamics of detuned quantum systems. The formalism also works for ensemble- 
averaged dynamics of stochastic systems. To illustrate the technique we give examples involving 
Raman processes, Bloch-Siegert shifts and Quantum Logic Gates. 
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I. INTRODUCTION 

In analyzing the dynamics of quantum mechanical sys- 
tems one is often confronted by Hamiltonian terms which 
are rapidly varying in time. And yet such off-resonant 
terms, which can require elaborate theoretical or numeri- 
cal handling, often end up having rather mundane effects. 
For example, consider a weak, far off resonant harmonic 
perturbation is applied to a two- level system (Fig.l). Its 
interaction picture Hamiltonian is 



H = 



Ml 



{|2)(l|cxp(-iAi) + |l)(2|exp(zAt)}, (1.1) 



where f2 is the Rabi frequency, which characterizes the 
strength of the interaction (and is assumed to be small 
compared with the detuning A). When we calculate 
the dynamics of such a system in the Bloch-vector rep- 
resentation [l[, we find a very rapid precession of the 
Bloch vector about the effective torque, which very nearly 
points along the z-axis. And yet, the effect of such a per- 
turbation is a slight change of the resonance frequency of 
the two levels (know as the "A.C. Stark shift" in quan- 
tum optics), described by an effective Hamiltonian of the 
form 

Ml 2 

^ = _ 4A {|2>< 2 I-|1><1|}- (1-2) 

This begs the question: is there a straightforward way 
to transform the rapidly oscillating Hamiltonian eq. (|1.1[) 
into a time averaged effective Hamiltonian eq. (|1.2p which 
accounts for the actual effect of the interaction? 



II. DERIVATION OF THE EFFECTIVE 
HAMILTONIAN 

In the interaction picture, the state of a quantum sys- 
tem (i)) evolves according to the following formula, 



where U (t, to) is the unitary time evolution operator 
which obeys the equation 



ih—U(t,t )=H I (t)U(t,to) 



(2.2) 



where Hj (t) is the interaction Hamiltonian. We have 
assumed the interaction picture here and throughout this 
paper. [l9| 

We are considering the coarse grained or time averaged 
dynamics of the system, defined formally by the following 
averaging procedure for some operator O (t) : 



0(t) 



f(t-t')d(t')dt', 



(2.3) 



where the function / (t) is real valued, has unit area, 
and is assumed to act as a low-pass filter (so that high 
frequency terms disappear from the average) . Time aver- 
ages defined in this way are well-known in describing the 
dynamics of quantum systems: the Rotating Wave Ap- 
proximation |lj is routinely used to discard various high 
frequency components in Hamiltonians of systems. Our 
time averaging is nothing more than a generalization of 
this approximation. Causality requires / (t) = if t < 0, 
although this is not a strict requirement of the theory. As 
is well known, the use of the rotating wave approximation 
leads to some non-causal predictions Q ; this will also be 
true of our generalization of the the R.W.A.: one can- 
not use these techniques in calculations where causal be- 
haviour is critical. Thus, while it may be strictly correct 
to use a causal averaging kernel for coarse-gained or time 

12) 



detuning A 



|1> 



\1>{t)) = U(t,to)\1>{to)) 



(2-1) 



FIG. 1: A two level system with a detuned interaction. 
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averaged dynamics, a non-causal function (e.g. a Gaus- 
sian function peaked at t = 0, and therefore non-zero for 
t < 0) may be easier to manipulate mathematically. In 
what follows we will avoid where possible making a spe- 
cific choice of / (t) , instead using generic properties of 
time averages. Using integration by parts, it can then be 
shown that the time derivative of the average is the av- 
erage of the time derivative, i.e. dO {t)/dt — dO (t) jdt. 

The effective Hamiltonian may now be defined 
straightforwardly by considering the time derivative of 
the time- averaged time-evolution operator, viz., 



ih^U(t,t )=H ef f(t)U(t,t ) (2.4) 

If we apply the averaging procedure discussed in 
eq. (|2.3p to eq. (|1.2[) and compare it with the eq. (|2.4p we 
obtain 



H eff (t)U(t,t )=H I (t)U(t,t ), 
which immediately leads to the formula 



(2.5) 



(2.6) 



Note that, although U (t, to) is unitary, the time aver- 
aged operator is in general not. Thus this expression 
Eq. (|2.6[) has the serious drawback that it is not Hermi- 
tian. The origin of this non-unitary evolution is easy to 
see: we have, in effect, divided our Hilbert space into a 
low- frequency part, and a high-frequency part; time av- 
eraging is equivalent to performing a partial trace over 
the high-frequency part. Performing a partial trace over 
one part of the Hilbert space is a classic paradigm for de- 
scribing non-unitary dynamics. One can carry out a pro- 
cedure entirely analogous to the derivation of the master 
equation (see for example Q) for this particular case, 
leading to the result that the effective Hamiltonian for 
the unitary part of the evolution is uniquely given by the 
Hermitian part of H. e ff (t), i.e. 



Heff (t) 



Heff (t)+H 



eff 



(2.7) 



We can obtain an expression for "Heff (t) using the 
standard expansion for the time-ordered evolution op- 
erator, U(t,to) = Texp -i y t Hi (*') dt'/h . Keeping 

terms to the second order of the Hamiltonian Hi (t) we 
obtain 



Heff (t) = Hj (t) + Hj (t) (t) - Hj (t) Ui (t) + 



where 



Ui(t) = 4 / H I {t')dt'. 



(2.8) 



(2.9) 



Since the time-averaging function f(t) was assumed to 
be real, it follows that the Hermitian adjoint and time 



averaging commute; using the fact that Hj (t) = Hj (t) 
and C7i (t) f = -Ui (t), we find that 

H eff (t)=W^j+l 



H(jb),U!(t) - H(t),U!(t) 

(2.10) 



III. HARMONIC TIME DEPENDENCE 

It is often the case that one deals with systems under- 
going a perturbation with an harmonic time dependence, 
at one or more frequency. Examples of this include an 
atomic or molecular system interacting with one or more 
laser beams; or an electron or nuclear spin interacting 
with an oscillating magnetic field; it is not difficult to 
think of other cases. In such systems the interaction 
Hamiltonian will have the following form: 

N 

Hi (t) =y^ h n exp (—iu n t) + h) n exp (iw n t), (3.1) 

71=1 

Where N is the total number different harmonic terms 
making up the interaction Hamiltonian. We assume 
(without loss of generality) that uj n > 0, and that 
^i < ^2 < • ■ ■ < wjv- We will assume that the range of 
frequencies present in the interaction is small, i.e. that 
max m ,„ \uj n - LO n \ = lu n - lui « lui. 

In this case, we can immediately calculate the operator 
Ui (*): 



Ui (t) = Vi (t) - V x (t ) 



where 



Vi(t) 



N 

E 

n=l 



1 



hi 



(3.2) 



(3.3) 



If the initial time t is in the distant past [20l ]. we may 

assume that V\ (to) = V\ (to), and immediately the term 
cancels from eq. (l2.10[) , implying that 



H eff (t) = Hj (t) + - 



(3-4) 

We will assume that the time average procedure acts as a 
low-pass filter. In other words, we are only going to con- 
cern ourselves with dynamic processes occurring at low 
frequencies; all processes oscillating at frequencies uj\ and 
above are assumed to average out to zero. Mathemati- 
cally, we assume that, for all values of m, n 

(3.5) 
(3.6) 
(3.7) 



exp (±iuj n t) 



exp (±i[cj„ 



exp (±i[uj„ - uj m ]t) = exp (±i[u) n - w m ]t) 
Using these assumptions, one immediately finds that 



Hi(t) 
Vi(t) 



(3.8) 
(3.9) 
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Substituting from eqs. (|3.7|) and (|3.9p into eq. (|3.4p one 
therefore finds, after a little bit of algebra, the following 
compact formula for the effective Hamiltonian: 



Heff (t) 



N 

E 



i 



h 



exp (lyuJn 



(3.10) 

where uj mn is the harmonic average of uj m and w„, viz., 



1 



2 V u) m 



1 



(3.11) 



This formula eq. (|3 . 1 0[) is the principal result of this pa- 
per. It was derived in a different manner by one of us 
a few years ago [l8[ , in a paper dealing principally with 
interactions between laser and trapped ions, and the ap- 
plicability of such systems as a possible implementation 
of a quantum computer. Since then various authors have 
employed the technique [f| @, 0], prompting us to un- 
dertake a more thorough investigation of the formula's 
applicability. To summarize, the assumptions that have 
led from eq. (|3.1|) to eq. (|3.10|) are (i) the interaction is 
sufficiently weak that we can we can terminate the series 
expansion of the time evolution operator after two terms; 
(ii) that the system is not responsive to high frequency 
stimulation, allowing us to make the time-averaging as- 
sumptions embodied in eq. (|3.5|) -eq. (|3.7|> : and (iii) the in- 
teraction has been occurring for a sufficiently long time 
that transient effects can be neglected. Note however 
that the theory is not a perturbation theory; if the effec- 
tive Hamiltonian is sufficiently simple that Schrodinger's 
equation can be solved , such a solution would embody 
all powers of the original interaction Hamiltonian. The 
solution of the propagator in these circumstances would 
be equivalent to a partial resumption of a Dyson series 
for the propagtor. 

Examinations of time-averaged dynamics have been 
presented by other authors. In particular, Avan et al. @ 
examined the problem using the Schrodinger picture, and 
obtained an expression involving sums over matrix ele- 
ments which is related to ours. More recently, the notion 
of an effective Hamiltonian for specific models of non- 
linear optical interactions has been presented by Klimov 
et al. b 9j. Furthermore, there is a considerable body of 
work in the nuclear magnetic resonance literature deal- 
ing with the notion of averaged Hamiltonians, which is a 
related concept (see for example |10j). However, to our 
knowledge, the compact formula eq. (|3.10p given here was 
not employed by any of these authors. We will now give 
some examples of its utility. 



IV. EXAMPLES 
A. The A.C. Stark Shift 

As discussed above, the A.C. Stark shift is an effective 
frequency shift experienced by a two level system due to 



the application of a far off resonance harmonic pertur- 
bation. As mentioned above the Hamiltonian describing 
this system is 

Hi (*) = ^{l 2 )( 1 l cx P(-* A ^) + |l)(2|exp(iAt)}. 

. (4.1) 

Immediately we can make the identification hi = 
(an/2) |2)(1| and ui = A. Applying eq. ([3~T0|) . and eval- 
uating the commutator, we quickly find that 



Heff = 



nn 2 
'1a 



{|2)(2|-|1><1|}. 



(4.2) 



Note that, if the detuning A were a negative quantity, 
we would have had to identify hi = {Ml/ 2) |1)(2|, since 
we have assumed u>\ is positive; carrying through the 
calculation results in a change of sign of the effective 
Hamiltonian, as one would expect. 

A related effect is the Bloch-Siegert shift [l|, which is 
a small shift in the resonance of a two level system which 
occurs due to the very high frequency components of the 
dipolar interaction, normally neglected under the Rotat- 
ing Wave Approximation. If one introduced the R.W.A. 
in a systematic way by employing some time averaging 
procedure, one could carry out an analysis analogous to 
that presented here, resulting in a small frequency shift 
term. One could also argue that the Lamb shift could 
be ascribed to an A.C. Stark shift caused by the vac- 
uum fluctuations of the electromagnetic field, although 
in this case one would have to deal with the on-resonant 
components in a different manner. 



B. Raman transitions 

As the next level of complexity, consider the case of 
a three level system interacting with two time-harmonic 
terms (e.g. fig. 2) 

In this case the interaction Hamiltonian is 



Hi(t) 



|3)(l|exp(-iAit) 



2 

+ ^|3)(2|exp(-*A 2 i)+J».a. (4.3) 



where h.a. stands for the Hermitian adjoint of the pre- 
ceding terms. Again, it is straightforward to make the 
identification hi = (KCli/2)\3)(l\, h 2 = (hn 2 /2)|3)(2|, 




FIG. 2: A three level system with two detuned interactions. 



4 



u>i = Ai and 0J2 — A2; substituting into eq. (l3.10|) and 
evaluating the commutators (four, in total), we find the 
effective Hamiltonian given by 



H, 



eff 



/,ri? {|3)(3|-|l)(l|}-^{|3}(3|-|2)<2|} 



4Ai 
4A 

-|2)(l|exp{-i(A 1 



4A 2 

l)(2|exp{t(Ai - A 2 )t} 
A 2 )t}). 



(4.4) 



The first two terms can be identified with A.C. Starks 
shifts associated with the two lasers; the second pair of 
terms represent transitions of population between the 
levels |1) and |2). Despite there being no direct inter- 
action between these two levels (indeed, for dipole inter- 
actions, if the transitions |1) — > |3) and |2) — > |3) are 
allowed, then |1) — ► |2) is forbidden), their population 
can be induced to oscillate coherently. These are the 
well-known Raman transitions, which play an important 
role in many areas of atomic and molecular interactions 
with light. 



C. Quantum A.C. Stark Shifts 

Consider the situation shown schematically in figure 
3. An ion is confined in an harmonic well, and is free 
to oscillate in the z-direction, and is interacting with a 
laser. Let us assume that the internal degrees of freedom 
of the ion can be modeled as a two level system, while the 
external degrees of freedom (i.e. its displacement from its 
equilibrium position in the trap) will be described by the 
standard quantum harmonic oscillator ladder of equally 
spaced quantum states. 

The Hamiltonian describing this interaction is as fol- 
lows 

Hi (t) = ^L\2){i\ e lk ^^- lAt + h.a.. (4.5) 

In this case A is the detuning between the laser and the 
internal resonance of the ion, k z is the component of the 




trapped ion 



harmonic trapping 
potential 

FIG. 3: A trapped ion, free to oscillate in the one direction, 
and interacting with a laser field. 



wavevector of the laser along the z-direction and the op- 
erator z{t) is given by 



k z z(t) = tj (ae~ tuJot + aV Woi ) 



(4.6) 



where 77, Lamb-Dicke parameter, is a dimensionless cou- 
pling constant, aJ and a are the creation and annihi- 
lation operators for the harmonic oscillator, and ujq is 
the trap frequency (see [ll[ for a more detailed discus- 
sion). Usually the dimensionless parameter 77 is small, 
hence we can make the "Lamb-Dicke" approximation, 
expik z z(t) « 1 + ik z z(t). This allows us to make the 
identification of three harmonic terms in this Hamilto- 



hi 
h 2 
h 3 



4m 



2>(l|ot 



2><i|a, 



(4.7) 
(4.8) 
(4.9) 



with u>i — A — loq, ll>2 — A and U3 
the effective Hamiltonian is 



A + usq. In this case 



H, 



eff 



V 



nn 2 

~ 4A~ 

x{|2)(2|-|1)(1|} 



(h + 1/2) , 



(4.10) 



where n = a' a is the number operator for the harmonic 
well. Once again, this effective Hamiltonian has the form 
of a frequency shift, but in this case the shift is dependent 
on the population of the the harmonic oscillator. Such 
'quantum' A.C. Stark have been studied before [12| : in- 
deed, since one may, using an appropriately detuned laser 
pulse, coherently map the internal degrees of freedom of 
an ion onto the external degrees of freedom, interactions 
of this kind have been proposed for use as quantum gates 

nam. 



D. Quantum gates 

Finally we consider the application of our formula to 
the derivation of quantum gate operations (this discus- 
sion first appeared in our earlier paper [18:] in which 
eq. p.lOp was first derived). Quantum gates between 
two-level qubits based on trapped ions and employing 
off-resonant interactions wereproposed by M0lmer and 
S0rensen a few years ago [HI, [l(| • These proposals have 
the distinct advantage that they perform well even with 
ions in any mixed state of the ions' collective oscillation 
modes, and were subsequently used to perform the de- 
terministic entanglement of four ions [17|]. A laser field 
with two spectral components, detuned equally to the red 
and to the blue of the ions' atomic resonance frequency is 
applied to a two ions in the trap. In the Lamb-Dicke ap- 
proximation, the interaction is described by the following 
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Hamiltonian: 

Hj(t) = ^—J(+) {1 + iri (ae-^ at ^-fite*" *)} cos(At) 

+h.a. 

= ^ J. - ^e^dt J, 
2 2 



In this equation — J x — i,J y is the collective raising 
operator for the ions (i.e. the sum |2)(1|^ + |2)(1|b of 
operators for the two ions A and B) and, as before, A 
is the detuning of the laser beam from the resonance 
frequency of the two-level system. Applying our formula 
eq. p.lOp for the effective Hamiltonian we find 



H, 



nn 2 v 2 

4(5 + u x ) 

any 

4(5 - lu x ) \5 + lu x 



J yd^Qj^ J y 

2lo x 



Titt 2 ri 2 
4(5 -u x ) 



Qi Jyj JyCt 



(4.12) 



This interaction is equivalent to a conditional quantum 
logic gate preformed between the two ions, and can be 
used to create multiparticle entangled states. 



V. CONCLUSION 

This paper has presented a simple derivation of a useful 
compact formula, eq. (|3.10[) . for the effective Hamiltonian 
for course-grained or time- averaged dynamics. We have 
found this formula to have been of considerable utility in 
a number of situations, and have presented some exam- 
ples here. We would welcome any comments from readers 
on other applications of this technique. A more detailed 
exposition, with a more rigorous analysis of the validity 
of the formula is in preparation [la ]. 
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